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Abstract 

We revisit the problem of magnetic field generation in accretion flows 
onto black holes owing to the excess radiation force on electrons. This 
excess force may arise from the Poynting- Robertson effect. Instead of a 
recent claim of the generation of dynamically important magnetic fields, 
we establish the validity of earlier results from 1977 which show only small 
magnetic fields are generated. The radiative force causes the magnetic 
field to initially grow linearly with time. However, this linear growth holds 
for only a restricted time interval which is of the order of the accretion time 
of the matter. The large magnetic fields recently found result from the fact 
that the linear growth is unrestricted. A model of the Poynting-Robertson 
magnetic field generation close to the horizon of a Schwarzschild black hole 
is solved exactly using General Relativity, and the field is also found to 
be dynamically insignificant. These weak magnetic fields may however be 
important as seed fields for dynamos. 

keywords accretion, accretion disks — galaxies: active — plasmas — magnetic 
fields — stars: magnetic fields — X-rays: stars 



1 Introduction 

The classical battery mechanism of magnetic field generation is connected with 
a noncoincidence of surfaces of constant pressure and constant density, where 
forces connected with pressure gradients become nonpotential or rotational. In 
this situation no static equilibrium in the gravitational field is possible. When 
considering separately the motion of electrons and ions, there is always a differ- 
ence in the velocities of electrons and ions which creates electric currents and 
an associated magnetic field. Self-induction is very important in the battery 
mechanism, determining the rate of increase of the magnetic field. 
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Along with ion and electron pressure gradients, a nonpotential force field 
may arise due to the radiation force which acts predominantly on the elec- 
trons. In a spherically symmetric star, the radiation force has a potential so 
that no magnetic field is generated: equilibrium in the two-fluid plasma results 
from a distribution of an electric charge and a static radial electric field. For 
geometrically-thin, optically-thick accretion disks, Bisnovatyi-Kogan and Blin- 
nikov (1977; hereafter BKB) showed that the radiation force above the disk 
has a nonpotential or rotational component. Under this condition, no electric 
charge distribution can give a static equilibrium. Instead, electric currents and 
a corresponding magnetic field are generated. The radiation forces above a thin 
disk gives rise to poloidal electrical current flow and a toroidal magnetic field. 

In accretion flows at very low mass accretion rates, an optically-thin, geo- 
metrically-thick accretion flow is possible (Shapiro, Lightman, & Eardley, 1976) 
where the ion temperature is close to the virial temperature. In the absence of 
a magnetic field, and neglecting relaxation processes between electrons and ions 
except for binary collisions these flows are referred to as advection dominated 
accretion flows (ADAF) (Ichimaru 1977; Narayan & Yi 1995). In the ADAF 
regime the radiative efficiency of accretion may be very low, ^10"^ times less 
then the standard value for a geometrically thin, optically thick accretion disk. 
Account of processes connected with the presence of a magnetic field increases 
the efficiency up to at least 1/3 of the standard value (Bisnovatyi-Kogan & 
Lovelace 1997, 2000). Nevertheless, the disk remains geometrically thick in the 
optically thin regime due to high ion temperature. 

Contopoulos and Kazanas (1998; hereafter CK) proposed that a cosmic bat- 
tery may operate in ADAF accretion flows owing to the Poynting-Robertson 
(PR) effect. The PR effect acts to generate a toroidal electrical current and 
poloidal magnetic field. The authors found that the magnetic field may be am- 
plified up to ^ lO^G in the vicinity of a black hole of stellar mass. Note that 
the PR mechanism of magnetic field generation is similar to the mechanism of 
BKB based on the nonpotential radiative force, where the magnetic field reached 
values of < lOG for a stellar mass black hole. In an optically thin disk both 
mechanisms act together leading to the generation of toroidal and poloidal com- 
ponents of the magnetic field. The influence of the PR effect on the dynamics of 
the surface layer of an accretion disk was treated by Mott and Lovelace (1999). 

Here, we analyze the difference in conclusions between CK and BKB. The 
radiative force initially causes the magnetic field to grow linearly with time. 
However, this linear growth holds for only a restricted time interval which is 
of the order of the accretion time of the matter. In CK the interval of linear 
growth is unrestricted. Even though we conclude that the magnetic field due 
to the radiation force is weak, it may have a role as a seed field for an a — w 
dynamo (sec for example Colgate and Li 2000). 

Section 2 treats the generation of toroidal field for the case of a thin disk, 
while §3 treats the generation of poloidal field in an ADAF flow. Section 4 
gives a General Relativistic treatment of a simplified model of PR magnetic 
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field generation in an accretion flow close to a Schwarzschild black hole. An 
appendix gives an explicit solution to the non-relativistic induction equation for 
the magnetic field generated by the PR effect in an ADAF flow. 



2 Radiatively induced current and toroidal 
magnetic field production in accretion disks 

Above a geometrically thin accretion disk around a black hole, the electrons are 
acted on by a nonpotential radiation force F l due to Thomson scattering. This 
was calculated by BKB, 

Fl - -i?COSe V0L = (i^LT?,, i^LS. 0) , (1) 

where a spherical coordinate system (i?, 9, <p) is used, and (f)L is the "pseudopo- 
tential" of the radiation force, which may bo expressed as 



ar f°° H{r)rdr 



c (i?4 + r4 + 2i?Vcos26l)V2 



(2) 



Here, the disk thickness is neglected and the cylindrical radius is r = RcosO. 
The function H{r) is the radiative flux emitted per unit area from one side of 
the disk. In the standard local accretion model, 

3 GMM ^ 

^W = ^^^- (3) 

where J = I — (rin/rY^^ , and Tin = 5rs = QGM/c^ is the inner radius of the 
disk for a non- rotating black hole of Schwarzschild radius rs- In the disk plane, 
9 = 7r/2, the radiative force is perpendicular to the disk, 

f,, = -^i^(.) = 5^^-^J, (4) 

where 

^ GMM ^ AncGMmp 

^ = , ^Edd = , (5) 

where ar is the Thomson cross section. Due to interaction of the radiation flux 

mainly with the electrons, the accretion disk becomes positively charged up to 
a value where the electrostatic attraction of the electrons balances the radiation 
force. The vertical component of the electrical fleld strength Eg in the disk 
plane is written as 



Eo{r) = -^Mr) = -E,-±- (^) V , (6) 
c|e| Lsdd ^ r J 
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where 

Thus, the surface charge-density of the disk is 

Eeir) 



27r 



The influence of this charge on the structure and stabihty of the accretion disk 
is neghgible. 

Both the gravitational and electrical forces have a potential, so that they can- 
not balance the nonpotential radiation force. Due to the radiation and electric 
forces, electrons move with respect to protons, which to a first approximation do 
not acquire the poloidal motion. Thus a poloidal electrical ciirrcnt is generated 
with an associated toroidal magnetic field. The finite disk thickness may create 
poloidal motion of all of the matter of the accretion disk, similar to meridional 
circulation in rotating stars (Kippcnhahn & Thomas 1982). The absence of this 
circulation occurs for a unique dependence of the rotational velocity over the 
disk thickness, Cl{z). 

To estimate the magnetic field strength, we write the electromotive force 
(EMF) as 

£ = ^ (j)FL X dl=^ J J dS xFl Egh , (8) 
where h is the half-thickness of the disk. Thus the stationary current-density is 

(Te£ (TeEgh 
Jst ^ ~ , (9) 

r r 

where Oe is the conductivity of the disk plasma. The stationary toroidal mag- 
netic field (BKB) is 

47r 47rcre Egh? 

B^o ^ — Jsth . (10) 

c c r 

In the radiation-dominated inner region of the standard a-disk model, h can be 
written as 

h = 3 J J rin , (11) 

(Shakura 1972; Shakura & Sunyaev 1973). Finally, we obtain the stationary 
value of toroidal magnetic field in the disk. 
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We next discuss the value of conductivity ae ■ 

Bisnovatyi-Kogan and Blinnikov (1977) considered different values of the 
plasma conductivity Ue, namely, the conductivity owing to binary collisions, 
o'couh and the effective conductivity Ce// derived by Vainshtein (1971), 

acoui^3xlO'T'/h'\ a,ff^^§^. (13) 

V -KCmO 

The magnetic Reynolds number in a turbulent plasma is defined as 

i?e™o = 5 , (14) 

where the turbulent velocity Vt in an a-disk is vt — acs, where Cs = \ppTp- 

In addition to the values of equation (13), we consider the conductivity in the 
presence of well-developed turbulence (Bisnovaty-Kogan & Ruzmaikin 1976), 

2 

^ ^ coul /-I -\ 

O-turb = -j T— = S ■ (15) 

47ran.Cs -Kcmo 

Estimations of the stationary field B^q of equation (12) in the radiation 
dominated inner region of the disk around a stellar mass black hole gives values 



(BKB) for two cases in ^ of B^q - 10"G for = acoui and B^q - lO^G for 
Of. = (Te//, with Re^Q ^ 3 x 10^°. For a turbulent conductivity ae = ^turb, we 
obtain iJ^g ~ 10 G for a = 0.1. 

The time-scale r™ for reaching the stationary field given by equation (12) is 
determined by the self-induction of the disk. This is equivalent to the "£ over 
time" of circuit with inductance C and resistance TZ. This time-scale is equal 

to 

Tm - 5 • (16) 

The crossing time of matter passing through the radiation dominated region of 
the disk is 

tc~- . (17) 

During the time tc there is linear growth of the magnetic field after which the 
matter falls into the black hole. Thus the stationary value of the large scale 
toroidal magnetic field is 



n ^ ( ^ \^ f" 
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B^^B^o— -3-£;o(-^) ['—) . (18) 

Tm Vr \LEddJ V r / 

For the case of a turbulent conductivity Ue — cFturb, the growth time-scale 
of the magnetic field, using equations ( |l5|) and ([l^), is equal to 



turb 



r 



T--«— . (19) 
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Taking into account that Vr — acs{L/LEdd){rin/r) < aCs and i*"''^ < tc, we 
find 

B, = « 3— So (j^) r-^Yj' . (20) 

aCs \LEddJ ^ r J 

In that 

s \LEddJ \ r / 

The strength of the stationary toroidal magnetic field produced by the battery 
effect in the radiation-dominated region of an accretion disk with the turbulent 
or higher conductivity from ( p^ ) is equal to 

At a distance r — 3ri„, we have 

2 /Mq 



This agrees with the findings of Bisnovatyi-Kogan and Blinnikov (1977). The 
corresponding magnetic energy-density is very much less than the energy-density 
associated with the turbulent motion in the disk pw^/2. 

3 Production of a poloidal magnetic field in op- 
tically thin accretion fiows by 
Poynting- Robert son eff'ect 

In optically thin accretion flows (Shapiro, Lightman, & Eardely 1976; Ichimaru 
1977; Narayan & Yi 1995), the radiation flux interacts with the inspiraling 
matter by the Poynting- Robertson (PR) effect (Robertson 1937). Analysis by 
Shakura (1972) showed that the PR effect was negligible for optically thick 
accretion disks. Contopoulos and Kazanas (1998) studied the PR effect as a 
mechanism for generating poloidal magnetic fleld in an optically thin accretion 
flow. They concluded that dynamically important magnetic field strengths could 
result from this effect. Here, we reconsider the PR effect for quasi-spherical 
advection dominated accretion flows (ADAF, Narayan & Yi 1995). 

The linear growth of the magnetic field due to the radiative force on the 
electrons found by CK is similar to that analyzed by Bisnovatyi-Kogan and 
Blinnikov (1977), but the PR effect implies an additional (small) numerator, 
{V(f,/c). Also, for a quasi-spherical accretion flow the characteristic scale is r 
instead of h, and the quasi-spherical luminosity is L/(47rr^) instead of H in 
(||). Then, using equations (p^), (p^, (p^, we obtain the rate of growth the 
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poloidal magnetic field due to the PR effect, whicli is equivalent to the expression 
obtained by CK, 



B.^B.,±^^(-I-)(^y('-^). (23) 
r,„ 3a \LEddJ \ r J \ r J 

Here, Eq is defined in equation (^, and in equation (|^). 

Now it is essential to take into account, that an element of matter with the 
induced magnetic field reaches the black hole in time tc ~ r/vr- (Damping 
of the magnetic field as it approaches the black hole horizon is discussed in 
Appendix A.) This means that the magnetic field grows only during the time 
tc- Consequently, the magnetic field reaches a stationary value 

3a LEdd V r / a LEdd M \ r J 

Taking into account that the luminosity is < IQ^^LEdd for optically thin accre- 
tion, and taking a = 0.1, we get a maximum value of the magnetic field created 
by the PR effect in an ADAF to a black hole, 

« 7 X 10-3 (^^) G . (25) 

This estimate of the field is about 10 orders of magnitude less that the value 
obtained by CK. The difference in the estimates results from the fact that CK 
assume that magnetic flux accumulates continuously near the black hole during a 
long time, reaching the equipartition with the kinetic energy. The accumulation 
actually occurs only during the time the plasma (which carries the field or 
current loops) takes to move inward to the black hole horizon (Bisnovatyi-Kogan 
and Ruzmaikin, 1976). The current loops created by the PR effect disappear as 
the matter approaches the horizon (see Appendix A). In the case of accretion 
onto a neutron star or a white dwarf, matter containing the current loops merges 
with the stellar matter, which is typically much more strongly magnetized. After 
merging, the matter becomes optically thick, the action of PR effect stops, 
penetration of matter into the magnetosphere of the star occurs, and interaction 
of the accretion flux with the stellar surface takes place. 



4 Magnetic field Generation in the vicinity of 
the Schwarzschild radius 

Here we consider the Poynting-Robertson magnetic field generation on the ac- 
cretion flow near a Schwarzschild black hole. The accreting matter is assumed 
to radially towards the black hole with velocity Vr- As before we consider the 
case of high conductivity matter where tm 3> tc- In a non-accreting flow, the 
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magnetic field can grow linearly with time as accepted by CK. As mentioned 
above, there is an important relativistic effect close to the black hole: Current 
loops in the accreting matter which approach the horizon of a black hole cannot 
produce a magnetic field visible by an external observer. This effect is related 
to the damping of the magnetic field in a collapsing star (Ginzburg & Ozernoi 
1964). 

The azimuthal force and the corresponding azimuthal electric field due to 
the Poynting-Robertson effect are 

Lax V^o . n p(ph) _ 1 pPR fr,fy\ 

^^^-wv^^' ^1' ' ^^^^ 

where 

v^o = A^I^, A^l, (27) 

and where the (ph) superscript indicates the physical value of the field com- 
ponent. Note that in a strictly spherical accretion flow there is no azimuthal 
EMF. However, in the approximate model considered here the infalling matter 
has a small rotation so that the PR force is small and these have radial inflow 
only slightly. 

We assume a Schwarzschild metric, 

ds^ = gooc^dt^ + gndr"^ + g2ide'^ + 5(33#^ , (28) 

where 

.oo=(l-^) , 5n = -(l-^)" , 

2 2 • 2 /I 

522 = -r , 533 = -r sm 6 , 

and Tg = 2GM/(?. The matter is free- falling in the radial direction with the 
nonzero components of a 4— velocity 



(l - ^) \ t^o = 1, (29) 



In any 4— dimensional time-space with metric gik (Latin indices takes values 
0, 1, 2, 3), the electric and magnetic i?" fields (Greek indices run over the 
values 1, 2, 3) in 3— space are defined through the antisymmetric electromag- 
netic field tensor Fj^ = —F^i, with zero diagonal components(Landau & Lifshitz 
1988), 

= , = Fo„ , (30) 
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and 

where Sap-y = e"^'' is the usual antisymmetric 3— dimensional symbol (£123 = 
1), and 7 is the determinant of the 3— dimensional metric tensor, obtained by 
splitting the metric tensor gik into space (70/5) and time (h) parts as 

laf} = -gap + gOaQop/gOO , h = .goO • (31) 



For the Schwarzschild metric (^8|), jap — —gap- The first pair of Maxwell 
equations can be written as 

Fik;l + Pli;k + Fkl-i — 

Note that the physical r and 6 components of the magnetic field in this reference 
frame are ^—guB^ and \/~g2iB^ respectively. Thus the dimensions of are 
length X B'^ . 

In a perfectly conducting medium moving with 4— velocity u', we have 
Fik'u^ = 0, which corresponds to a vanishing electric field in the comoving 
frame. This gives 

Fa = -^f~g Eap^iv^ lc)B'' ■ (33) 

In the presence of an externally imposed electric field i?^^*, the electrical field 
Ea is 

Fa = -V^ eapyiv^/c)B^ ~ El^' . (34) 
The 3— velocities are given by 

- ^"^^^ dx" = cdi, (35) 



,j2\ -1/2 

72 



where 



c 



1 / 2\-i/2 



Substituting the expression ( ^ ) into equation (^2|) , gives the following equa- 
tion for the magnetic field i?" 

' ^r-^iB^v-_B.,,)^^,.,.^_Xc. (37) 



dt 8x13''' " Q^p 
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For the Schwarzschild metric ( p8[ ) {x^ , , , x^) — {ct,r,9,(j)), it follows from 
(pgh, @, @ that = K',0,0) and 



r Mr 



(39) 



where j c? — Tgjr. The value of /i, the determinant g of the 4— metric tensor 
(7ifc , and the determinant 7 of the metric tensor 7q,^ in the Schwarzschild metric 
are 



h=\- 



g = r sm 



-1/2 



sin t 



(40) 



With = {B^B^Q) and = (0,0, E^^^*), and with all quantities 

independent on the azimuthal angle 0, equations (|37|) and (pS) give 



di^/^v'-B^) 
Ft 



.diy^v'-B^^ 



Or 



dEf 
dr 



iVlB") = 



(41) 

(42) 
(43) 



9r'^ ' '89 

Equations (^) and (^) with known right-hand-sides are solved using the char- 
acteristic method. 

The integrals of the characteristic equations can be written as 



t 



dr 



Jrr 



dr 



-gv'^B'' + c 



06 v-Vh 

Q^ext 



C2 



dr 



-dr = C. 



(44) 

(45) 
(46) 



The constants Ci are determined by the initial conditions. For the present 
problem these are 



= B" 



0, r = 7'o at i = 



which implies Ci — 0. In the general case the constants C2, C3 are determined 
by the initial values of B^ , B^ , which should satisfy zero divergence condition 
([43I). We may in general take Ci = 0, fixing the reference frame r = ro at i = 0. 
With account of equations (^9|) and (^o|) , equations ( ^5|) and (^^ can be written 
as 



d 



r^/^y/r^ sin dO 



(47) 
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(48) 



Integration of equation (|4J) with Ci = gives a relation between r and 
ro in the form 



ct 
r„ 



3/2 



+ 2x1''^ + In 



V^- 1 

1 

1 



(49) 



+ 1 



where x — r/vg, and x^ = r^/rg (Bisnovatyi-Kogan and Ruzmaikin, 1974). 
Consider now the Poynting- Robertson EMF, E^^*{r,6). First we show that 



£^5^' must tend to zero as 



at least as fast as {r 



or faster in order to 



avoid singularity at . This means that in the comoving reference system with 
metric 



(50) 



there is no singularity at the black hole horizon. The connection between 
Schwarzschild and comoving coordinates (r, p) (the angle coordinates 9 and 
6 are the same) is 



CT 



ct 



ct 



2y/x + In 



1 



+ 1_ 



In 



(51) 



We can now connect the magnetic field in a comoving system B" = {BP, B^ ,0) 
with the field in the Schwarszschild system in terms of Schwarszschild variables 
(r, t) as 

BP = , (52) 



B' 



B' 



r„ 



El 



r r{r — rg) saiO 



It follows from equations ([47|), ([iq), and ( ^4 ) that there is no singularity in 



the comoving frame if £"5^* tends to zero as [r 



or faster as r 



The 



metric tensor in this system ( |50| ) is regular on the horizon, so with finite B°' 
all 4— invariants (for example, FikF'^^) are also regular there. In fact, we can 
obtain the dependence of E'r^^ from equation (p9), taking into account that 



the covariant component E'^ in equations (|41|) and (^2|) is connected with the 
physical component E^F^^ from equation (p6|) as 



hE^ -rsmOE^ 



(53) 
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The luminosity L seen by a distant observer viewing collapsing matter with 
constant comoving luminosity Lq is 



L = io 1 - ^ 



(54) 



(Zeldovich & Novikov 1971). Thus we have from equations (p6[), p7|), (p3|), and 
5l 



where 



L> = 



sm 



LoOtAs/GM 
47rc2|e| 



(55) 



(56) 



Equation ( |55|) is of course simplified but the dependence allows an estimate to 
be made of the magnetic field generation by the PR effect close to a black hole. 
It is necessary that -E^'^'* vanish sufficiently rapidly on the horizon in order to 
avoid a physical singularity, but the exact dependence is not important. 
Substituting equation (|5^) into (Ell) and (B^) gives 



2D ^/F^ 



In 



ro 



2\rJ r 3\roJ 2 \ro J 



^0. 



B" = 



D 



TgT^r - r. 



3/2 



1 



1-^ 



2 V^o 
1 

4 

sin( 



cos 9 



(57) 



' g 



(58) 



There are several limiting cases where expressions for B^ and B^ can be written 
in a simpler form. 

1. Non-Relativistic, Newtonian Regime: Here, r, ro ^ rg, and from 
equation (H9) we have 



3 ct 
2a;3/2 ~a 



Xq = X [ I + 

and from equations (^) and (|8|), 

B^ = l^ln(l 



2/3 



B" = 



Dsinl 



1- 1 



3 ct 

2.T'V2 ^ 



3 ci 



2x3/2 



(59) 

(60) 
(61) 
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We see here that for large t the physical value rB^ tends to a finite limit of 
the order of equation (|2^) , while B^ grows but only logarithmically. During the 
accretion time to a massive black hole this logarithm does not exceed ^ 25. Thus 
the magnetic field is larger by a factor ~ 25 than the estimations of equations 
( p4[ ) and (|25|). Still, the magnetic field is enormously less than the value found 
by CK. 

2. Vicinity of the Gravitational Radius: Here, (a; — 1) <C 1 and we have 
from equation ( ^9| ) 



1 = 4- 



1 



+ 1) 



ct 



CXP + -Xq 



3/2 



2.T, 



1/2 



(62) 



For matter initially situated in the vicinity of the horizon the Lagrangian 
coordinate a;o ^ 1 ^ 1- In this case we have 



In 



Xq-I 



ct 



From equation (57) 



B'' = 



D cost 



2r. 



5/2 



(a; -1)3/2 [exp(4ct/rg) - 1] 



From equation (58), 



B*" 



Dsm0 

7/2 



(63) 



(64) 



(65) 



These relations are valid also at large t while {x — l)e^*/'"« <C 1. 

For matter with an intermediate Lagrangian coordinate, {xq 
vicinity of the horizon r 



B'' = 2D cos 



g is reached only at very large t, so that 



5/2 



B" 



11 1 
L>sin6' {xo-iy 



2xq 



_3_ 

Xo 



7/2 



1 



11/2 



1) - 1, the 



(66) 



(67) 



We see that in the vicinity of the horizon (x — 1) <C 1 for matter with an inter- 
mediate Xq, {xq — 1) ~ 1, the component B^ in the Schwarzschild coordinates di- 
verges. However, this is a natural coordinate singularity which means only that 
a Schwarzschild observer cannot exist physically in this region. Consequently, 
the magnetic field in Schwarzschild coordinates has observable consequences 
only in the region x ^ 1, r ^ rg. 

Matter at a very large Lagrangian radius reaches the vicinity of the gravita- 
tional radius at very late times. Consider a case with 



Xq — I :$> X, 



1 < 1. 



(68) 



13 



Here, for very large ct/rg '3> \ ln(x — 1)| we have a;o — (3ci/2rg)^/'^ and 



AD cos 9 , ,1/9 / 3ct \ , . 

5^ = — T7^(^-l)'/'lnU^l (69) 



If we assume, in addition that xo >> (x — 1) wc obtain 



B'^^{x-l)y\ (70) 



Dsine 

'9 

From a comparison of the asymptotic relations (^) - (^9|), we see that for 
r approaching r^, the Schwarzschild component grows exponentially with 
time very close to the horizon while remaining zero at the horizon. For larger 
r it grows logarithmically which is the dependence in the Newtonian domain. 
The Schwarzschild component has the mentioned singularity on the horizon, 
which is unobservable because a physically realizable observer cannot measure 
it. However, at any fixed value of r close enough to r^, the temporal behavior of 
the B^ can be obtained from equations (|65|), (|67|), and ([70|). In this region B^ 
starts to grow exponentially with time, but at long times it tends to a constant 
value. This component rapidly decreases with increasing r. 

It is of interest to determine the magnetic field "seen" by an observer lo- 
cated far from the black hole who measures the field near the gravitational 
radius by means of the cyclotron radiation coming from this region. The ob- 
server determines the field strength by measuring the cyclotron frequency of the 
radiation. The frequency of the emitted radiation is determined by the comov- 
ing field strength, that is, by B^ and & from ( p^ evaluated in the vicinity 
of Tg. We now obtain estimates of these fields. It follows from equation ( |62[ ) 
that matter with Lagrangian coordinate xq approaches the horizon as t — > oo. 
The coordinate ccq parametrizes the horizon points and the radial coordinate x 
parametrizes points of the initial hypersphere t = Q. From equations (|5^), (p7|), 
(|6^), and we obtain the comoving fields in this region. 

It is easy to verify, that for any value of rp , < tq < oo , we have 

< ^y^lnxo. (73) 

For accretion to a stellar mass or massive black hole the logarithmic factor does 
not exceed ~ 25, and the value of B^ remains of the same order of magnitude 
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as the Schwarzschild component in the Newtonian region, equation (|6^) at 
large times, formally extrapolated to r ^ r^. The poloidal component , given 
by equation (|7^), is equal to zero at ro — rg, and ro — oo, and has a maximum 
at r = llrg/3. Thus 

, ^fl, AZ? sinfl ,_ . 

B' < \B'Ua. = , (74) 

where A = (3^8^11"^^)^/^ w 0.04. Therefore, the possible values of B^ are less 
than the Newtonian value of the Schwarzschild component B^ (|6|) at large time 
near rg. 

The proper cyclotron frequency loq of radiation emitted at (r, t) , but mea- 
sured with respect to the comoving time is 



LOo — 

mc 



"^'{BPf + r^B^f 



1/2 

(75) 



where m is the electron rest mass. Using equations ( [73| ) and (Q), we obtain an 
estimate for the upper limit on this frequency. 



D\e\^/2 I , 

[uJojniax ~ gyj y In Xo + — ■ (76) 

mcrg * * 

The frequency measured by a distant observer lu is is related to ujq as 



--oV^f^^^, (T7) 
1 — V cos V'/c 



where v/c = \/rg/r and is the angle between (which is in the nega- 
tive r-direction) and the direction of the photon trajectory in Schwarzschild 
coordinates. For a radially emitted photon {ip = tt) we have from equation 
([77[),tJ = 0^0(1 ^ -^/r^/r) near the horizon, and for the tangential direction 
(jp = 7r/2) we have w = wo(l ~ ^g/^)- Due to the upper limit on a distant 
observer does not see the light emitted very close to the horizon, because of the 
very large red shift. Thus for a distant observer the relativistic region close to 
the horizon is unobservable. 

Thus we conclude that the magnetic field produced by the PR effect close 
to the horizon of a black hole can be safely estimated using the Newtonian 
approximation given by equations (|60| ) and ( |6l| ) . The estimated magnetic fields 
are dynamically insignificant. 

5 Conclusion 

We have reconsidered the battery effect in accretion flows due to the non- 
potential nature of the radiation force on the electrons. We considered cases of a 
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geometrically-thin, optically-thick disk where a toroidal magnetic field is gener- 
ated, and a geometrically-thick, optically-thin ADAF type accretion flow where 
a poloidal magnetic field is generated due to the Poynting- Robertson effect. For 
a stellar mass black hole the generated toroidal field is estimated to be < 10 G, 
while the poloidal field in an ADAF flow is < 0.01 G. The fields vary inversely 
with the black hole mass. In both cases the flelds are dynamically insignificant. 
The very large fields obtained by CK resulted from assuming unrestricted linear 
growth of the magnetic field. The field grows only during the accretion time. 
A General Relativstic treatment of the Poynting- Robertson generated magnetic 
field close to the horizon of a black hole shows that that the field magnitude 
may be larger by a factor < 25 than the values obtained with a non-relativistic 
treatment. Even though the magnetic field due to the radiation force is weak, 
it may have a role as a seed field for an a — w dynamo (see for example Colgate 
and Li 2000). 
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Appendix A: 

Equation for Poloidal Magnetic Field Due to PR 
Effect 

Here, we treat in more detail the influence of the azimuthal radiation force F^^ 
which is rotational^ and it cannot be balanced by any axisymmetric electrostatic 
field. The radiation is mainly from the central region of the flow so that the 
radiation flux-density \s S L/{AtiR^), with L the accretion luminosity, and R 
the distance from the origin. The PR radiation force on an electron is = 

—SoTVcji/c^, where ctt is the Thompson cross section and is the azimuthal 
velocity of the accreting matter. Including the radiation force. Ohm's for the 
plasma is 

J = CTe(E™ + E + vxB/c) , (78) 

where is the electrical conductivity, and E-^^ = cf) S<7TVfi,/{\e\c^) is the 
Poynting-Robertson electric field. 

Combining Faraday's and Ampere's laws and equation (^Sl) gives 



^^^ + v.VvI/ = cri?|'« + ,;eA'^*, (79) 

where rje — /('^■ncFe) is the magnetic diffusivity, and ^ — rA^ is the flux 
function and is the toroidal component of the vector potential. Also, A* = 
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q2 jQ^2 _ (\jr)dldr + /dz^ in cylindrical coordinates and A* — /dR^ + 
[(1 — y?)/ R^]d'^ /dy? in spherical coordinates where ^ — cosO. Note that 
Br = —{l/r)d'i>/dz and Bz — {\/r)d'^ /dr in cylindrical coordinates, while 
Br = {R'^sinOy^d^ /d9, and Be = -{Rsmey^d"^ /dR in spherical coordi- 
nates. Taking v^, = {GM/Ry/^g{0), with g{Tr/2) = 1 and g{0 ^ 0, tt) = 0, we 
find 



mpc^g L _ 2gEoJ^ /r^^5/2 

63/2|e|Li5rfdi?5/2 ^ 63/2 L^^^ { ji 



{^y\ (80) 



where Eq[M) is given by equation (|7|). Equation (79) is equivalent to equation 
(8) of CK. 

For an ADAF flow, the poloidal velocity is vr — —a^{GM/R)^^'^ where 
a is the Shakura-Sunyaev parameter and ^ < 1 is a constant (Narayan & Yi 
1995). Following Contopoulos and Kazanas we write rj^ = VRIvrI, where V, the 
magnetic Prandtl number, is the ratio of magnetic diffusivity to viscosity. Mea- 
suring R in units of ri„ and t in units of to = r^n /{a£, V GM) = (ri„ /c) / (a^) , 
equation ( |8C| ) becomes 

Kgi9)sm9 1 ^ ^{9'^ , I - f^' d^^\ 

dt~ i?3/2 i?2 5^2/' ^ ' 

where K = rl^Eo{L/LEdd)/{?'aS,). 

The time-scale of the linear growth of ^ is Tm = to/V. For a turbulent 
magnetic diffusivity where V = 0{1), this time-scale is quite short, ~ to ~ 
7.3 X lO"VM/M0)(O.l/a)(l/^). Therefore, the physically relevant solution to 
equation (|8l| ) is the stationary one where the Poynting-Roberston term cx if is 
balanced by diffusion. This gives Kg{6)sm9 = —V{1 — ii^)^^^ jdp? so that 
is independent of R. For example, for giQ) = sin(^?), Br ~ (K/V) cos9/ R^ or 



B 



0.6 cos 6* L M /rj„\2 



PB. 



a £,V LEdd Mq 



and Bg = 0. This estimate agrees with equation (24). Equation ( p^ ) corresponds 
to a radially outward field in the northern hemisphere and a radially inward field 
in the southern hemisphere. The polarity of the field agrees with the PR drag 
on the electrons in the —(f) giving a ring current in the +(f) direction, while the 
simple nature of the field results from the approximation that 5* w L/{AitR'^). 
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